Ferromagnetic phase transition and Bose-Einstein condensation in spinor Bose gases 
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Phase transitions in spinor Bose gases with ferromagnetic (FM) couplings are studied via mean- 
field theory. We show that an infinitesimal value of the coupling can induce a FM phase transition 
at a finite temperature always above the critical temperature of Bose-Einstein condensation. This 
contrasts sharply with the case of Fermi gases, in which the Stoner coupling I s can not lead to a 
FM phase transition unless it is larger than a threshold value lo- The FM coupling also increases 
the critical temperatures of both the ferromagnetic transition and the Bose-Einstein condensation. 
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The magnetism of Fermi (electron) gases has long been 
a research topic in solid state physics. Although many 
open questions remain, the magnetic properties of Fermi 
gases have been well understood 0. The Fermi surface 
plays an important role in determining their magnetism. 
For example, a magnetization density M in Fermi gases 
increases the band energy by splitting the Fermi surfaces 
for spin- up and spin-down particles. As a result, Fermi 
gases mainly behave as Pauli paramagnets in the absence 
of an exchange interaction. If an effective ferromagnetic 
(FM) exchange I s is present, electron gases can exhibit 
ferromagnetism. Within the framework of the Stoner the- 
ory, I s results in a negative molecular field energy when 
M is finite. If I s > Iq, the Stoner threshold, the value 
of the molecular field energy becomes larger than that of 
the increase of the band energy induced by M, and then 
a FM ground state is energetically favored. 

The magnetism of Bose gases has been less studied. 
But since the realization of Bose-Einstein condensation 
(BEC) in ultracold atomic gases[2|, more and more at- 
tention has been attracted to this topic, because the con- 
stituent atoms, such as 87 Rb, 23 Na, and 7 Li usually have 
(hyperfine) spin degrees of freedom and thus a magnetic 
moment m. m can be large in atoms such as 52 Cr[3j, for 
which m = 6fJ-B, where \ib is the Bohr magneton. Since 
1998, atomic gases can be confined and cooled in purely 
optical traps in which their spin degrees of freedom re- 
main active, and therefore investigating their magnetic 
properties becomes experimentally possible |J|. 

Ferromagnetism in spinor bosons without any spin- 
dependent interactions has already been theoretically 
studied. As opposed to the case of fermions, the ground 
states of spinor bosons are degenerate and a ferromag- 
netic state is among the ground states 5, 6] . Furthermore, 
the spinor Bose gas is rather apt to be magnetized by an 
external magnetic field even at a finite temperature T, as 
long as T < T c , the BEC critical temperature^ H R. lib] . 

In ultracold atomic g Heisenberg-like exchange 

interaction is usually present, arising from spin-flip scat- 
tering processes. For example, the effective two-body 



interaction in spin-1 atoms was given bv|llj|. 
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where ra(r) = Y, a V4( r )^( r )' S a (r) = 
E CTCT < ^iM-S^/i/v-lr) with S a (a = x,y,z) being 
3x3 spin matrices, Vv( r ) is a field annihilation operator 
for an atom in the hyperfine state \F — I, a) at point 
r, and a = 1,0,-1. The normal ordering removes the 
spin-spin interaction of a particle with itself. In the 
special case of 87 Rb, the exchange is ferromagnetic 
(c s <0)E3 El- 
Ferromagnetic couplings can also be generated by "su- 
perexchange" processes in correlated spinor boson sys- 
tems described by the boson Hubbard model. Recently, 
this model was used to describe Bose gases moving in 
a periodic optical lattice[l4|, and a lot of work has been 
done to discuss its properties^!!- It has been shown that 
the spinor boson Hubbard model can be reduced to the 
ferromagnetic Heisenberg model, provided the on-site re- 
pulsion is large enough 6] , 
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(2) 



where Ih is the Heisenberg exchange constant. 

The presence of FM couplings enriches the properties 
of spinor bosons significantly. One can naturally imagine 
that a FM phase transition could be induced, at a tem- 
perature plausibly depending on the energy scale of the 
coupling. Thus, an interesting question arises as to how 
the FM transition interplays with the BEC. We focus on 
this question in this paper. Without loss of generality, we 
study a homogeneous spinor Bose gas before considering 
trapped gases. Our theoretical results indicate that the 
FM phase transition takes place at a finite temperature 
Tp always above T c , regardless of the size of the coupling. 
This contrasts sharply with the case of fermions. Figure 
1 shows schematically the relation between Tp and I s 
for both Bose and Fermi gases. More recently, Isoshima 
et aZ-E3| and Huang et a/.0] studied BEC in trapped 
F = 1 spinor bosons with FM couplings, but the FM 
transition was not considered. 
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FIG. 1: Schematic relations between transition temperatures 
and FM couplings I s . T F and T' F represent the FM transition 
temperature for Bose and Fermi gases respectively. T c and To 
denote the BEC critical temperature for spinor bosons with 
and without couplings respectively. Iq is the Stoner threshold. 



The model - Following the Stoner theory for fermion 
gases 0, a molecular field H m = I S M is introduced to 
describe the effective FM coupling between bosons, and 
the energy shift arising from the molecular field is 



H m S 1 



-I,MSf 



(3) 



The molecular field contains all the spin-dependent inter- 
actions. It can be derived from the exchange interactions 
via the mean-field approximation. In order to avoid the 
spin-spin interaction of a particle with itself, we treat the 
particles as being on some kind of "lattice" , 



(4) 



where Z is the effective "coordination number" , which for 
a gas is an irrelevant dimensionless parameter of order 
unity, and J\f is the number of sites. Therefore I s — 
ZIh for the Hamiltonian in Eq. (2). M, = (S^) serves 
as the ferromagnetic order parameter. We investigate a 
spinor Bose gas with hyperfme spin-F and set the boson 
magnetic moment to unity. It is convenient to choose 
(Si) = (0,0, Mi) and ~M l = (Sf) = Ea^iM wherc 
annihilates a boson with spin quantum number a at 
site i. We take Mi = M for a homogeneous boson gas. 

Our main purpose is to show how the FM coupling 
affects the properties of spinor bosons. So we neglect 
spin-independent interactions for simplicity. Then the 
effective Hamiltonian for the grand canonical ensemble 
reads 



H — N/x = ^[fk - M - f C#m + H e )]nk a , 



(5) 



where ek = H 2 k 2 /2m* is the kinetic energy of free parti- 
cles with mass m*, \i is the chemical potential, and H e 
is the external magnetic field. Since the Hamiltonian is 



diagonal, we may calculate the grand thermodynamical 
potential 

n = -hnZ = -4mTre-' 3 ^- Ar ^, (6) 
P P 

where Z = Tr{exp[— (3(H — Nfj,)]} is the partition func- 
tion, and = l/(fcsT). The density of particles is 



where V is the volume of the system, N is the total 
number of particles, and n = N/V. M is determined 
self-consistently by 
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(8) 



The external magnetic field H e is set to zero in the fol- 
lowing calculations. 

We now consider only F = 1 bosons. Eqs. (7) and 
(8) lead to the basic equations determining the phase 
diagram of the spin-1 boson system, 



1 = no + t ^ 
M = n^+ti 



/a(o) + /i(o + 6) + /a(o + 26) ,(9a) 



h(a)-fs(a + 2b) 



(9b) 



where a = -(/x + H m )/(k B T), b = H m /(k B T), the 
reduced temperature and coupling are given by t — 
k B Tm*/{2Trh 2 n 2 ^) and / = I s n l ^m* / '(27rft 2 ), respec- 
tively, uq is the condensate density, no — n^/n is the 
condensate fraction, M = M jn is the normalized magne- 
tization, and fsifl) is the polylogarithm function defined 
as 18] 
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P 



(10) 



We note that / s (0) = C( s )? the Riemann zeta function. 

Eq. (9a) is the standard formula for the density of 
bosons, by which the reduced BEC critical temperature 
t c can be determined [l9j. a > for t > t c and a — > 
as t — > t c from above. Hq > and a = for t < t c . 
In deriving Eq. (9b), we assume that only the spin-1 
component of the bosons can condense. This assumption 
is discussed in the following. 

Preliminary analyses - Assuming a FM phase transi- 
tion is induced by I below the reduced transition temper- 
ature it is of special interest to determine the relation 
between and /. We first suppose I is very large so that 
tp > t c - Provided that the FM transition is continuous, 
i.e., 6^0 with t — * tp, Eqs. (9) become 



1 = 3t F fs(ap), 
1 = 2l4h(a F ). 



(11a) 
(lib) 
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FIG. 2: Reduced temperature t vs. reduced FM coupling 
/ phase diagram of spin-1 Bose gases. The paramagnetic 
normal phase (PM), ferromagnetic normal phase (FM), and 
coexisting phase of ferromagnetism and Bose condensation 
(FM+BEC) are indicated. The solid (t c ) and dotted (t F ) 
curves are results from Eqs. (15) and (11), respectively, and 
diamonds (t*) and circles (t F ) are obtained from Eqs. (9). 
The insets show the differences between t* and t c , and be- 
tween t* F and tp, respectively, versus /. 




FIG. 3: Plots of the normalized magnetization M c at reduced 
BEC critical temperature t c vs. reduced FM coupling /. The 
solid curve denotes M c at t c obtained from Eqs. (15), and the 
diamonds represent M c at t* obtained from Eqs. (9). Inset: 
plots of tp — t c (inverted triangles), and tp —t c (solid curve). 
The dotted curve is a guide to the eye. 



at the BEC critical temperature, the self-consistent Eqs. 
(9) reduce to 



where ap = a(tp). Equations (11) define a relation be- 
tween tp and I. For a given tp, I is given by 



3/a (a F ) 



2U 2 (a F ) 

J is a monotically decreasing function of ap. As a 
/|(a F ) -» C(3/2) « 2.612, and fi(a F ) » ^Ja~p~ 
for small values of I and oj?we have 

ap Antol 2 , 



(12) 

-+0, 
. So 

(13) 



where t = l/[3C(3/2)] 2 / 3 is the reduced BEC critical 
temperature for the Bose gas with 7 = 0. Since aF is 
always larger than zero as long as I is finite, Eq. (13) 
implies that an infinitesimal FM coupling can induce a 
FM phase transition at a finite temperature above the 
BEC critical temperature, because of is finite only for 
t > t c . 

We now calculate the asymptotic behavior of the rela- 
tion between tp and I for small couplings: / <C 1. From 
Eqs. (11) we find[i3 



t_p 

to 



— « 1 + 8TTtll. 



(14) 



This equation shows that the FM transition temperature 
increases with the FM coupling. 

By decreasing t further, the Bose gas will undergo 
BEC. Due to the molecular field, the spin-1 bosons have 
the lowest free energy. So the condensate contains only 
the spin-1 component and the other two boson spin com- 
ponents remain in their excited states. Since no = 



1 = ti 



JL — +2 



C(2)+/|(M+/|(26c) 



(15a) 
(15b) 



where b c — b(t c ). These equations define the relation 
between t c and I. Since the FM phase transition occurs 
above t c , the system has a finite magnetization at t c . 
When I< 1 and b c is very small, Eq. (15a) leads to 



tc 

to 



4(l + V2)V^b~c 



(16) 



Equation (16) is similar to the expression of the BEC 
critical temperature in an external magnetic field0, but 
here b c w &irt[)I 2 is not the external field, but is a 
self-consistently determined quantity proportional to the 
spontaneous magnetization. Substituting b c into Eq. 
(16), we find 



5< 2 



2)TTtnI. 



(17) 



This equation indicates that the FM coupling also in- 
creases the BEC critical temperature. 

Numerical results- Comparing Eqs. (14) and (17), 
we see that tp is slightly smaller than t c for small FM 
couplings. This is unphysical and might bc attributed to 
the assumption that the FM transition is continuous. To 
investigate this point, we solved the self-consistent Eqs. 
(9) numerically. The results obtained indicate that both 
the FM transition and the BEC become discontinuous if 
the FM coupling is sufficiently small. 

Fig. 2 shows the phase diagram of spin-1 Bose gases. 
We first discuss the Bose-Einstein condensation. For 
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FIG. 4: Order parameters of FM phase (M) and Bose con- 
densed phase (no) versus the reduced temperature t for the 
reduced FM couplings / = 0.5, 0.4, 0.3, from top to bottom. 
The diamonds mark the points where BEC takes place. At tp 
(circle), M drops discontinuously from M F = M(t F ) to zero. 



so close to t* that it is very difficult to solve Eqs. (9) for 
the FM normal phase. 

It is worth noting that the FM transition is contin- 
uous for large couplings. This point is consistent with 
the Weiss molecular field theory for classical particles, 
in which the FM transition is continuous Q, because for 
large couplings, the FM transition occurs at a relatively 
high temperature, when the Bose statistics reduces to 
Boltzmann statistics. 

In conclusion, we studied the ferromagnetic phase tran- 
sition and Bose-Einstein condensation in spinor Bose 
gases with ferromagnetic couplings via mean-field theory. 
We showed that the coupling, regardless of its magni- 
tude, induces a ferromagnetic phase transition at a tem- 
perature always above the critical temperature of Bose- 
Einstein condensation. Moreover, the ferromagnetic cou- 
pling also increases the critical temperatures of both 
phase transitions. 

R.A.K acknowledges partial support from the Max- 
Planck-Institut fur Chemische Physik fester Stoffe. 



I > I c w 0.2, t c = t* c , as shown in the left inset, and 
both a and no are exactly determined to be zero at t c , 
while M has a finite value: M(t c ) = M c . In this case, 
the results from Eqs. (15) agree perfectly with those ob- 
tained from Eqs. (9). With increasing /, t c (t*) tends to 
the upper limit value f /[C(3/2)] 2 / 3 which is the reduced 
BEC critical temperature of scalar bosons. Nevertheless, 
when I < I c , the BEC becomes a discontinuous phase 
transition: no has a small but nonzero value at t*. In 
this case, t* is slightly larger than t c , as the left inset 
shows. The fact that the Bose gas has a finite magneti- 
zation at t c (t*) implies that a FM transition has taken 
place at a higher temperature as stated above. Fig. 3 
shows the normalized magnetization M c at t c (t*) versus 
the reduced FM coupling I. 

Fig. 4 plots the order parameters of both phases. 
The normalized magnetization M drops quickly with t 
increasing above t c {t* c ), and tends to zero at the re- 
duced FM transition temperature t* F . t* F = tp when 
J > I F « 0.35, as shown in the right inset of Fig. 2. 
The smaller / is, the more quickly the magnetization 
drops. When I < Ip, M does not disappear continu- 
ously but drops from a finite value Mp — M(t* F ) ^ 
to zero abruptly at t* F , with t* F > tp as shown in the 
right inset of Fig. 2. This signals that the FM transition 
becomes first order. Mp « O.f f for / = 0.3. In this case 
Eqs. (f f) no longer hold. We determine t* F for the cou- 
pling down to / = 0.22 and the results show that t* F > t c , 
although tp < t c when / < 0.24, as shown in the inset 
of Fig. 3. Hence the FM transition occurs at a higher 
temperature than the BEC. However, for small /, t* F is 
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